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^ . Abstract 

^ , Nanoscale optoelectronics and molecular-electronics systems operate with current injection and 



< 



nonequilibrium tunneling, phenomena that challenge consistent descriptions of the steady-state 

transport. The current affects the electron-density variation and hence the inter- and intra- 

^ . molecular bonding which in turn determines the transport magnitude. The standard approach 

^ . 

for efficient characterization of steady-state tunneling combines ground-state density functional 

C^ ■ theory (DFT) calculations (of an effective scattering potential) with a Landauer-type formalism 

and ignores all actual many-body scattering. The standard method also lacks a formal variational 

C \ basis. This paper formulates a Lippmann-Schwinger collision density functional theory (LSC- 

^ ' DFT) for tunneling transport with full electron-electron interactions. Quantum-kinetic (Dyson) 

CN ' equations are used for an exact reformulation that expresses the variational noninteracting and 

•^ , interacting many-body scattering T-matrices in terms of universal density functionals. The many- 

in 



body Lippmann-Schwinger (LS) variational principle defines an implicit equation for the exact 



t^^ ■ nonequilibrium density. 

o 

OO 

O ■ PACS numbers: 73.40.Gk,71.15.-mJ2.10.-d 
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I. INTRODUCTION 

The function of heterostructure-based optoelectronics and of future molecular electron- 
ics^ rests on current-injection and tunneling transport that causes genuinely nonequilibrium 
conditions. The systems are nanoscale and the performance is sensitive to the atomic con- 
figuration. The technologies define difficult theory problems of calculating nonequilibrium 
tunneling in the presence of electron-electron/boson scattering. A quantum cascade laser- 
(QCL) produces optical transitions in repeated stages with current injection in resonant- 
tunneling structures. The electron-electron interaction is ubiquitous, causes strong many- 
body scattering out of equilibrium,- and directly affects the QCL operation. This is be- 
cause electrons which exit one QCL stage must be recycled for current injection in subse- 
quent stages and because the optical activity depends on details of the energy distribution 
of injected electrons. Similarly, a current-driven molecular-electronics switch*^ positions a 
nanoscale molecule between leads and uses charge-transfer processes to adjust the intra- 
and inter-molecular bonding and morphology. Understanding details of the current- induced 
changes in the interacting electron distribution is also here important because the molec- 
ular morphology (and nature of bonds) determines the magnitude of the nonequilibrium 
tunneling. Deriving a quantum-physical description that provides ab initio, predictive and 
self-consistent accounts of nonequilibrium tunneling with full electron-electron interaction 
is very desirable. 

For equilibrium systems it is possible to rely on the traditional, ground-state density func- 
tional theory-i^ (DFT) to provide a materials-/system-specific account. The ground-state 
DFT is formulated in the canonical ensemble, conserves the total number of electrons, and 
takes as input only the potential defined by the nuclei. The approach uses predefined approx- 
imations for a universal functional that expresses contributions of the electron-electron inter- 
action to the total, interacting ground-state energy for a system in equilibrium. Traditional 
implementations of ground-state DFT calculations, using the local density approximatiort^ 
and/or the generalized-gradient approximation,- provide a rich description of the bonding 
inside materials with a dense electron-density distribution and within molecules. Recent ex- 
tensions^ii^iii of the density functionals to include accounts of dispersive interactions allow 
descriptions of vdW-DF bonding and organization in sparse materials^^ and in (as well as of) 
macromolecules.~ The ground-state DFT calculations define an ab initio level (as opposed 



to a model level) of details in the description of niaterials-/niolecular-physics properties and 
are extremely useful for they allow a transferable description of specific systems. It is this 
ab initio level of detail that we seek also for the open, nonequilibrium tunneling systems. 
Unfortunately, the openness and transport invalidate the assumptions of traditional DFT 
which rests in a ground- state-total- energy variation principle that applies in equilibrium. 

There is exciting recent progress towards ab initio calculations of steady-state interacting 
tunneling even if consistent quantum-physical calculations of open nonequilibrium systems 
are challenging. The nonequilibrium tunneling problem is difficult because the finite ap- 
plied bias makes phase space available for actual electron scattering processes (which are 
normally suppressed in equilibrium). It is also difficult because it is essential (but hard) 
to ensure conservation rules^^ in accounts of the nonequilibrium and interacting electron 
dynamics. The Landauer-type formulations,— characterizing tunneling transport in terms 
of ballistic propagation of individual particles moving in an effective potential, constitute a 
simple approximation. They can formally be derived^»i^»i^ for the linear-response regime or 
in the absence of many-body scattering using nonequilibrium Green functions^i^U (NEGF) 
within the so-called partition scheme introduced by Caroli et al,"^^ but they only have limited 
applicability. One solution strategy for ab initio calculations of nonequilibrium interacting 
tunneling involves use of quantum-kinetic equations--- for conserving calculations and there 
exists a number of formal results,— '^^'^Si'^i^^i^^ extending and correcting the Landauer-type 
formulation. Several explicit calculations for nonequilibrium tunneling with electron-electron 
interactions focus on correlated-electron systems like the Kondo problem and use both di- 



agrammatic approaches (for example, Refs. ll7lJ26lJ27l ) and exact reformulations^^i^^ of the 



quantum-kinetic equations. Tunneling through molecular systems has been investigated 
on an ab initio level diagrammatically starting from either Hartree-Fock eigenstates of the 
molecular system^^ or by applying a conserving GW formulation.— A second solution strat- 
egy for interacting steady-state tunneling invokes time-dependent density functional theory^ 
(TD-DFT) either in combination with a master-equation approach for a finite, closed-loop 
systent^i^ or with the NEGF formulation in an infinite, open, partition-free formulation 
with an explicit time-evolution in the applied bias.— i^ Traditional TD-DFT contains no 
implicit dephasing and describes coherent evolution of the (interacting) many-body wave- 
function for finite systems. Explicit and implicit dephasing mechanisms must be carefully 
analyzed^^i^ in these applications to steady-state tunneling. Both of the TD-DFT-based 



methods allow calculations of interacting tunneling transport that are in principle exact, 
and the TD-DFT basis should make it simpler to achieve an ab initio level of details. All of 
the above-mentioned solution strategies for ab initio nonequilibrium calculations are com- 
putationally intensive. 

The standard method for efficient calculations of steady-state tunneling with an ab ini- 
tio level of detail combines ground-state DFT with a Landauer-type formulation^ com- 
puting tunneling transport as an independent-particle transmission problem. The result- 
ing ballistic-transport DFT (BT-DFT) approach was introduced by Lang2^ and uses the 
ground-state DFT exchange-correlation energy to define an effective potential for scattering 
of independent particles. The BT-DFT represents a natural (but ad hoc) extention of the 
Poisson-Landauer-type transport solvers.-^ It is a meaningful approximation because elec- 
tron conservation is automatic in the ballistic single-particle description. There exist efficient 
implementations^^i^E and BT-DFT calculations provide valuable theory characterizations of 
tunneling systems. However, the BT-DFT approach also constitutes an uncontrolled approx- 
imation. Use of a ground-state DFT formulation must be discussed even in linear response.— 
Use of the Landauer-type formulation, ignoring all actual many-body scattering events, is 
problematic for fully nonequilibrium conditions. Analysis^ of the frequency-resolved current 
density shows that nonequilibrium electron-phonon scattering causes distribution changes 
that are inconsistent with Landauer-type descriptions (corresponding problems for nonequi- 
librium electron-electron scattering can be inferred from Ref. y). The use of a Landauer-type 
formulation prevents BT-DFT from consistent calculations of current-induced changes in the 
electron distribution and hence of effects that are important for nanoscale optoelectronics 
and molecular electronics. The fact that BT-DFT lacks a formal variational basis motivates 
a continued search for efficient ab initio calculations of nonequilibrium tunneling with full 
electron-electron interaction. 

This paper formulates a Lippmann-Schwinger collision density functional theory (LSC- 
DFT) for steady-state, nonequilibrium tunneling systems treated in the Born-Oppenheimer 
approximation. The LSC-DFT is based on formal collision theory^ for the interacting many- 
body problent^i^ and allows an exact discussion, resting on the Lippmann-Schwinger (LS) 
variational principle.— i^i^ The LSC-DFT is expressed through universal density functionals 
that characterize the variational form of the noninteracting and interacting many-body T- 
matrices. The LSC-DFT provides a formal solution in terms of an implicit equation for the 



exact electron density. The formal LSC-DFT solution constitutes a natural starting point 
in a search for rigorous formulations of single-particle schemes based on the LS variational 
principle. 

This paper is organized as follows. Section II defines the partition scheme and the Hamil- 
tonian. It also discusses the general (nonequilibrium and interacting) tunneling problem as 
a complex many-body collision problem. Section III presents the formal density-functional 
basis for the theory while section IV develops the LS T-matrix functional description. Fi- 
nally, section V contains a summary and outlook while appendix A provides details of the 
uniqueness-of-density proof for the time-dependent interacting tunneling problem. 

II. PARTITION SCHEME, HAMILTONIAN, AND COLLISION THEORY 

It is convenient to utilize the partition framework of Caroli et al, Ref. [2lj, but retain the 
full level of details of atomistic, many-body calculations everywhere close to the tunneling 
region r ~ at all physically relevant times t ~ in the collision problem. For simplicity, the 
tunneling structure only comprises a left (right) lead r < Zi{> zr) plus a central tunneling 
region 'C in 2;^ < r < zr. Atomic units will be used throughout and the full kinetic energy 
is written 

5^ y rfr' y dv ^+(r) (-\{V5{v - v')Y\ 4(r) ^KL + Kn + Kc + 5K (1) 

where Kl,r,c has a basis strictly confined to separate components and where 

5^= E E /^r^.^(r) (-\{V6{v-z,)Y\ur). (2) 

i=L/R s -^ ^ ^ 

At time t -^ —00 the partition scheme assumes that each of the disconnected subsections 
of the partitioned system Hd = Yli=L c r-^^^^ ^^ equilibrium at (generally) different chemical 
potentials Hl/c/r- The operators Nl/c/r describe the electron count in each of the three 
subcomponents. Initially, at t ^ —00, the system is described by a static potential fo(r) 
and corresponding operator 

Vo= [drvo{r)n{r) (3) 



where n(r) = ^^, ^+(r)^s(r) denotes the electron-density operator. I assume, for simplicity, 
that fo(r) reduces to a uniform background potential (pL/R (with a value set by the average 



electron concentration) far in the leads. The choice of initial Hamiltonian, 

i=L/R/C 

ensures an automatic charge neutrality at t ^ — oo (and at z ^ ±00 even at finite t). The 
equilibrium distribution at t ^ —00 is specified by a Gibbs energy weighting Hd — fiL^L — 
IJ'rNr — Hc^c which is independent of the value of the applied bias 0bias = f^L — fJ'R and 
which exclusively depends on the initial electron concentration in the leads (and on the 
initial electron occupation of the central island 'C'). 

The LSC-DFT further assumes an adiabatic turning on of the tunneling, the electron- 
electron interaction W, and of the static electron-scattering potential fsc(r) which includes 
the effects of the applied bias and of the set of atomic potentials. One can also allow for an 
addition of a time- dependent potential 4>g{r, t > to) that describes a possible gate operation 
starting at some finite time to- The time-dependent collision potential is 

Wcoi(r, t) = [t;sc(r) + 0g(r, t)] exp{r]t) (5) 

and the collision problem is formally described by the Hamiltonian 

H{t) = Ha + Hi{t) (6) 

Hiit) = [6K + W]exp{r]t) + 6V{t) (7) 

Vcoi{t) = 5V{t) = f dvVcoi{r,t)h{r) (8) 

where the factor exp{r]t),ri — > 0^ expresses the adiabatic turning on. The collision term 
that drives the dynamics is Hi{t). One may, without lack of generality, assume that the 
collision potential fcoi(r, t) also contains an implicit multiplicative factor that restricts the 
basis to a finite, but very large region (much larger than zi, < z < z-^\ given the choice for 
VQ{\r\ — > 00) this assumption simply amounts to treating the remote part of the leads as 
jellium. 

The expectation value of the electron density is defined^^ 

„(M) = (ftW>W.Iii|!g^ (9) 

nH{r\t) = f/(t,-oo)+n(r)f/(t, -00) (10) 



by establishing the initial (equilibrium) density matrix f){t -^ — oo) = po and by formally 
solving for the ensuing (nonequilibrium) dynamics described by^ 

i?M^=Hmt,-ooY (11) 

This applies to general quantum-statistical problems but it is important to provide con- 
sistent (conserving) approximations to the time-evolution f/(t, — cxo) of the systems. Such 
approximations are notoriously difficult to obtain for interacting systems out of equilibrium. 
In this paper, I use the LS many-body collision theory^i^i^ and describe the interact- 
ing tunneling problem as a complex collision problent^i^ in which (an ensemble of) fully- 
interacting many-body electronic states of the leads encounter and scatter off a potential in 
some central tunneling region. For any initial many-body state |$g) (eigenstate of Hd) one 
may formally obtain the many-body collision stated 

|vl/^,+(t))=?7(t,-oo)|$e) (12) 

from a direct application of the temporal-evolution operator ( TTTl) . The case of a purely 
static scattering potential is then described as an elastic many-body collision but the finite 
bias still causes actual electron-electron scattering^ that invalidates assumptions of ballistic 
transport. An inelastic collision event arises when the ensemble of many-body states scat- 
ters off a corresponding time-dependent collision potential.—'^ An effective time-dependent 
potential may also arise in the description of tunneling in the presence of a quantized boson 
field,— "^"^"^ as is relevant for further investigations of optoelectronic devices, for example, 
lasers under typical operating conditions. 

The many-body collision problem (for tunneling) is from the outset specified by the choice 
of partitioning, that is, (i) a specific choice of Hamiltonian Hd with an interrupted kinetic 
energy K — 6K and (ii) the lead electron concentration, and (iii) the choice of an initial state 
|$o) {oT ensemble of initial states, po) that also formally depends on the initial distribution 
in the central region 'C 

III. DENSITY FUNCTIONAL THEORY OF INTERACTING TUNNELING 

For a time-dependent collision potential it is relatively simple to generalize the TD- 
DFT analysis^*^ to the many-body collision theory of nonequilibrium interacting tunneling, 



appendix A. This demonstrates that the time- dependent density is a unique functional of the 
time-dependent coUision potential. The time-dependent collision density functional theory 
contains the steady-state formulation of LSC-DFT as a limiting case under some conditions. 

A. Density functional theory of time-dependent tunneling 

Applying a time-dependent gate voltage to a nonequilibrium tunneling system probes a 
response^*^ that reflects electron correlations. It is interesting in itself to develop a formal 
density functional theory basis for such time-dependent many-body collision problems. 

For a given bias and partitioning of the general time-dependent collision problem I analyze 
the mapping M : fcoi(r, t) — > nir, t) and argue 

01. The electron density n{r,t) is a unique functional of the collision potential Vcoi{T^,t). 

For the collision problem the goal is to determine the nonequilibrium time-evolution of a 
single state (or a grand-canonical ensemble) of an infinite open system under a collision term 
Hi{t) which includes an adiabatic turning on of both tunneling SK, the electron-electron 
interaction W, and a time-dependent potential 6V{t) (with a finite basis). In contrast, 
TD-DFT considers a single state which from the outset is described hj K + W and which 
evolves under an external potential V{t) and demonstrates uniqueness of the time- dependent 
density when the system is finite^ and/or when the external potential V{t) has a finite 
basis.— However, it is straightforward to generalize the reductio ad absurdum argument 
from TD-DFT to the present many-body collision formulation of interacting nonequilibrium 
tunneling described in a partition scheme. This is because the difference of collision terms, 
Hi(t)—H[{t), remains exclusively specified by the difference in collision potentials, fcoi(i"5 ^) — 
t>^jj[(r,t) (and this difference may be assumed to have a finite basis). 

Formally the observation 01 only establishes the uniqueness of the density variation and 
history for a given bias and for a given choice of partition scheme (that is, choice of 2;l/r 
and choice of the initial occupation in the central region 'C'). A different partition scheme 
produces a different unique mapping. A/" (as well as corresponding T-matrix functionals for 
the scattering behavior). However, the choice of partition scheme must become irrelevant for 
very large z-^ — zl and basis of the collision term 6V{t). This is argued on physical grounds 
for tunneling systems that lack singular responses: since dephasing eventually decouples the 
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time and spatial correlations^i the solution n(r, t) must eventually become insensitive to 
memory effects as well as details of the potential far in the leads.— i^ 

B. Density functional theory of steady-state tunneling 

For static collision problems the four- dimensional variational space of a general density 
history n(r, t) naturally becomes overcomplete in its definition of the scattering potential 
fsc(r). The mapping A/"^^ : n{r,t) — > Vco\{T^,t) shows that a general density history causes 
potential variations 



6v = — - — 6n{r,t) 



-1 



OTi 

Vr2(5r + -;:—St 
ot 



(13) 



5n Sn 

that are incompatible with the steady-state transport assumption. 

For time-independent collision problems that have a steady-state solution, n(r,t) = 
n(r, 0), I argue 

02. The time-independent scattering potential fsc(r) is uniquely determined by the steady- 
state density n{r) = n{r,0). 

The adiabatic turning on of the static scattering potential fsc(i") can be viewed as a limiting 
case of a time-dependent tunneling problem with the restricted variation: 

(^t-coi = Vfsc exp(?7t) 6r + r] v^c exp{r]t) 6t — > Vv^c St. (14) 

The unique mapping A/" : Vco\{r,t) — > n{r,t) identifies the corresponding set of relevant 

density variations. If the steady-state tunneling problem is characterized by non-divergent 

values of 6Af/6v we have^ 

dn SAf / \ ^ f.,^s 

-^ = ^ ^^sc exp{r]t) — > 0. (15) 

This is the condition that specifies the steady-state scattering solutions, Sn{r) = 6n{r, t = 0). 
The observation 02 permits formulation of universal density functionals that character- 
ize variational expressions for the noninteracting and interacting T-matrices in LSC-DFT. 
Formally these universal functionals also depend on the assumptions that are build into the 
partition scheme. For some tunneling problems it may be important to retain a functional 
dependence on the initial island occupation fie- Here I assume that the solution of the 
steady-state LSC-DFT problem is a functional only of the density (for a given applied bias 
and choice of ^l and fiji) . 



IV. A LIPPMANN-SCHWINGER T-MATRIX FUNCTIONAL DESCRIPTION 

Formulation of a single-particle scheme with independent dynamics of fictitious particles 
is important for an accurate and efficient evaluation of the electron density in interacting 
nonequilibrium tunneling problems. Prerequisites for such a formulation are (a) correspond- 
ing variational expressions of a many-body quantity evaluated both in the interacting and 
noninteracting cases, and (b) universal density functionals to characterize those variational 
physical quantities. The LSC-DFT uses the many-body LS variation principles for the inter- 
acting and noninteracting many-body T-matrices as well as exact reformulations to satisfy 
those necessary conditions. The LS variational properties also permit the LSC-DFT to 
specify an exact implicit equation for the nonequilibrium electron density. 

A. The Lippmann-Schwinger variation principle 

In their seminal paper^^ Lippmann and Schwinger identified a set of variational properties 
for the collision problems. The original theory ignores the self-energy shifts associated with 
the adiabatic tuning on of the collision term Hi (t) but the work was soon after supplemented 
and regularized by additional analysis.— i^ The LS variational principle^ applies for any 
combination of collision states and it should be straightforward to generalize the following 
also to finite-temperature tunneling problems. Like in traditional equilibrium DFT, however, 
the focus will likely remain on zero-temperature properties. It is natural to build a functional 
that refiects the evolution of the ground-state |$o) of the original disconnected system. 

The LSC-DFT provides an exact variational evaluation of the ground-state-to-ground- 
state transition matrix element Ti/[n], with usual definition Tv[n] = (<l'o|-f^i|^o,+ ['^])- This 
matrix element is a functional of the tunneling electron density because the scattering state 
|\E'o,+ M) is specified by the choice of external potential. Through the optical theorem,— 
this T-matrix element characterizes the total rate of tunneling (a charge transfer process) 
arising in the presence of full electron-electron interaction. 



The LS variational principle is expressed using the notation of Ref. 
the discussion in Ref. 



45 



42 



I generally follow 



and introduce G^ = {Eq — Hd ± irf) ^ as the retarded (+) and 



advanced (— ) Green function operator of the original disconnected Hamiltonian Hd while 

|^o,±N) = |$o) + Gli7i|^o,±N) (16) 

10 



identifies the many-body (collision) state that evolves forward/backwards in time under the 
collision term Hi. These states are functionals exclusively of the density n{r) (or n{r,t) in 
the wider problems beyond the present scope). I further introduce 

E.[n,V] = (M/o,-N|i/i|$o), (17) 

E4n,V] = ($o|i/i|^o,+ M), (18) 

T[n,V] = {^0,-HHi - HiGiHil^oAn]) (19) 

which, like the compensated form 

Tv [n]=E_[n,V]+E4n,V]-T[n,V], (20) 

are functionals of the density n but also contain an explicit dependence on ^ = Vq + 6V 
through the collision term Hi. All four functionals represent a correct evaluation of the T- 
matrix behavior when evaluated at the correct density n (the density n that results under the 



collision term Hi). For T[n, V] this fol 
T+ = l-HiGi = l- GiHi, see Refs. 



ows by a simple application of the Dyson equation, 



43 



45. 



The key observation for the LSC-DFT formulation is that the extremum identified by the 
variational condition 

6Tv[n] _ 5Ty 5^^,^ ^ 5Ty ^vl>o,+ _ ^ 

identifies the electron density n that solves the collision problem Hi. This follows from the 
(many-body) LS variational principle^i^ because the derivative (5Ty/(5\l/o,-(+) is directly 



45 



proportional to the many-body LS equation for scattering states |^o,+(-))) Ref. 

The noninteraction collision problem, defined by H^ = 6K + 6V, has a corresponding 
density functional description. It has a different unique mapping between the density and the 
potential 6V and different scattering states I^I/oi) ^^^ ^^i^ mapping defines other (related) 
functionals 

E'_[n,V] ^ {^l_[n]\H'i\%), (22) 

El[n,V] = ($o|ff?|<^M), (23) 

T°[n,V^] ^ {^l4n]\H', - H'.GiH'MAn]) (24) 

T^[n] = E'_[n, V] + El[n, V] - T^[n, V]. (25) 
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The extreinuin, identified by the variational condition 

^T°[n]=0, (26) 

identifies the density that solves the noninteracting problem H^ = SK + SV. 

B. Universality of T-matrix functionals in LSC-DFT 

To obtain a description given in terms of universal functionals it is necessary to identify 
the partial contributions that arise from the kinetic-energy addition, the electron-electron 
interaction, and the potential scattering and to find a method to evaluate the difficult many- 
electron effects once and for all. This is possible by formal manipulation using the Dyson 
equation and by use of the LS equation itself. By construction, the formal manipulation 
does not invalidate the underlying variational character of the LS formulation (l20l) . 

I first introduce scattering states and Green function operators for the set of partial 
collision problems defined at 6V = 0. I use |xo±) ^^^ \Xo,±) to denote the collision states 
for the 6V = noninteracting and interacting connected problems, given by Hd + SK and 
Hd + 5K + W, respectively. Also, G\ = [Eq - {Ha + 6K) ± tr]]~^ and G± = [Eo - {H^ + 
6K + W) ± iri]~^ identify the noninteracting and interacting Green function operator at 
6V = 0. These noninteracting and interacting collision problems contain an implicit choice 
of potential Vq while the general noninteracting and interacting collision problem is described 
hj V = Vq + 6V = Koi- Neither xo,± ^oi Xo ± are therefore functionals of n (and the same 
applies for all Green functions in use). 

For the matrix element S^[n] = ($o|('^-^ + '^^)l^o+M) I use a simple formal manipula- 
tion of the 'ket' state 

KM) = \<^o) + G'ii5K + 5V)K4n]) 

= \xl+)+GlSVK4n]), (27) 

as can be verified by a direct application of the LS equation.^' The resulting separation 

El[n,V] = {^o\SK\xl^) + {xl-\SVKM) (28) 

can, of course, be repeated for a separation also of E^ [n] and of the corresponding interacting 
expression 'E±[n, V], for example, 

E4n,V] = {%\{SK + W)\xo,+) + (xo,-|<5^|^o,+ W). (29) 

12 



This reformulation can also be derived, Eqs. (XIX. 9) and (XIX. 120) of Ref. |42|, by applying 
the Green theorem on the weighted overlap between the two collision states. The Green 
theorem plays a central role in NEGF calculations for the open tunneling systems and enters, 
for example, in the quantum-kinetic based derivation of resonant-tunneling rate equations.— 
For the noninteracting matrix element T°[n, V] = (\1'q _[?2]|((5-ft' -|- 5V){1 — G'^{5K + 
5V)}\^^ j^[n\) (and for T[n, V^]) it is necessary to first expand the 'bra' and 'kef collision 
states by the LS equation. I collect terms involving either a (xo-l oi' a (^o-MI 'bra' state 
and either a |xo+) or a l^o+M) 'kef state separately and I use the underlying quantum- 
kinetic (Dyson) equation for simplification. Taking one of the cross terms as an example, 
one obtains 

{xl^5K{Gl\V,] - G%bKGlm - G%W\^lM) = 0. (30) 

Repeated applications of the Dyson equation completes the separation 

-^'{n.y] = ixVmii - Gi6K}\xi^) + «_Mimi - G°+mK+N)- (31) 

A corresponding expansion applies, of course, also for the interacting matrix element 

T[n,V] = {xo,-\iSK + W){l-G'ii6K + W)}\xo,+) + {^o,-Hm'^-G+6V}\^oAn])- (32) 

The electron-electron interaction effects on the many-body scattering problem can now 
be expressed in universal functionals. A set of complex constants 

a_ = {xo,-\m + W)\<l>o) (33) 

a+ = {%\{6K + W)\xo,+) (34) 

b = {xo,-\{SK + W)[l-Gi{SK + W)]\xo,+) (35) 

(along with corresponding definitions a^, 6° for the noninteracting case) characterizes the 
dynamics in the absence of the collision potential (at 6V = 0). More importantly, a set of 
collision-state matrix elements 

A4n]ir) ^ (vI/o,-M|n(r)|xo,+) (36) 

A+[n]{r) ^ (xo,-|^(r)|^o,+ N) (37) 

Bi[n]{r) = (^o,-M|ri(r)|^o,+ M) (38) 

i32H(r,r') = (^o,-NI^(r)G+n(r')|^o,+ M) (39) 

13 



represents universal density functionals that determine the many-body dynamics when the 
coUision potential is included in the presence of full electron-electron interaction (while 
corresponding universal functionals ^^[n] and 'B°2M characterize the full collision problem 
at W = 0). The variational form of the interacting and noninteracting T-matrices can thus 
be reformulated 

Tv[n] = a. + a+-b+ f drv^dr) /CyH(r) (40) 

T^[n] = at + al-b' + j drv^r) /C°.[n](r) (41) 

JCP[n]{r) = ^f H(r)+^L°^[n](r)-Sf^M(r)+ f dr' Bf\n]{r,r')v,,{r'). (42) 

In essence, calculation of a set of universal functionals (for relevant choices of chemical 
potentials /Ul/r) permits a simple general evaluation of the interacting T-matrix for arbitrary 
scattering potentials fsc(r)- 

The LSC-DFT description also permits a succinct formulation of the interaction effects 
on the T-matrix functional derivatives 

6{Ty[n] - T^[n]) _ 6{]Cv[n] - ICUn]) ^ ^, . JAICv[n] 

The interaction effect is expressed as a complex function of r and it is entirely specified by 
suitable approximations to universal density functionals. 

C. Variational solution of the interacting collision problem 

Separating out the noninteracting dynamics (for which we can seek highly accurate char- 
acterizations) and the interaction effect ( l43l) defines a formal LSC-DFT solution 

5n(v) 5n(v) ^^ 5n(v) 

which constitutes an exact but implicit equation for the nonequilibrium electron density. 

The formal LSC-DFT solution (jHj) serves as a natural starting point for search for a 
single-particle scheme for calculations of the density in specific nonequilibrium tunneling 
systems. The possibility is exciting, for a rigorous single-particle scheme would permit 
efficient and exact calculations of (5Ty[n] and would ensure automatic current conservation 
in ab initio calculations specified by universal functionals. The single-particle LS equation 
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certainly determines the noninteracting many-body dynamics described by 5Ty[n]/6n. It 
is not a priori clear that the interaction term fH3l) represents an additional effect caused 
by some effective single-particle scattering and it is not a priori clear that a single-particle 
scheme exists for the LSC-DFT. However, the present results show that the LSC-DFT 
satisfies necessary conditions and motivate a search for rigorous single-particle formulations. 

V. SUMMARY AND OUTLOOK 

This paper formulates a Lippmann-Schwinger collision density functional theory (LSC- 
DFT) for nonequilibrium interacting steady-state tunneling. The theory rests on the 
Lippmann-Schwinger variational principle for the interacting and noninteracting many-body 
T-matrices and includes exact reformulations that express the variational T-matrix forms 
through universal density functionals. The variational property of the LSC-DFT specifies 
an exact implicit equation (jHj) for the electron density. The LSC-DFT furthermore fulfills 
necessary conditions for a possible formulation of a rigorous single-particle scheme. The 
present results motivate a future study (using the formal LSC-DFT solution (l44l) and the 
LS variational properties of single-particle scattering) to explore conditions on the dynamics 
and to test if a rigorous single-particle scheme can be defined in LSC-DFT. 

Of course, any implementation of a LSC-DFT method must also rely on successful for- 
mulation of a good approximation for the universal functionals ^±[n] and -61,2 [^] that char- 
acterize the complex many-electron collision behavior. The formulation of TD-DFT-based 
ab initio calculations^^i^^ facilitates a program to explore the T-matrix behavior for a range 
of scattering potentials and thereby deconvolute approximations for the universal function- 
als .4.±[r2] and i3i^2['^]- The partition-scheme method for time propagation of the tunneling 
many-body wavefunction^^ may allow a direct extraction of T-matrices and simplify the task. 
It is natural to first seek functionals that have a local-density flavor in the parameterization 
of Eqs. fl36|) - fl39|) but it is also possible that more complex functional forms must be explored. 
Exact solutions of nonequilibrium correlated-electron model systems^Si present possibilities 
for further refining parameterizations of approximations for the universal functionals. 
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APPENDIX A: UNIQUENESS OF DENSITY IN COLLISION THEORY 

I argue uniqueness of the time-dependent density 01 for the complex LS collision problem 
of open tunneling systems in the partition scheme as a relatively straightforward generaliza- 
tion of the TD-DFT analysis for finite systems^ and for infinite systems with a restricted 
basis of the single-particle potential.— The grand-canonical foundation, the use of the par- 
tition scheme, and the basis in quantum-kinetic equations makes the analysis of the many- 
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and the argument is included 



body collision problem slightly different from that of Ref. 
here for completeness. 

For a specific partition with given initial configuration and initial density matrix operator 
pit — > — oo) = po, I consider the time evolution 

^^-^ = [H{t\m], (Ai) 

with formal solution given by the many-body evolution operator 

pit) = U{t, -oo)^poU{t, -oo). (A2) 

I compare two similar systems given by Hi{t) and H[{t) for which the /c'th time-derivative of 
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The current-density 



the collision potential begins to differ at some time tj, Eq. (3) of Ref. 
operator 

j(r) = {2t)-' 5^[V^+(r)]^,(r) - V^+(r)[VV^,(r)] (A3) 

s 

constitutes a sensitive probe of system differences at times immediately thereafter, t = tf. 
The density matrices p{t) and p'{t) that correspond to Hi{t) and Hi{t) must, of course, 
agree at tj. Use of ( lAip permits the evaluation 

2-Aj(r,t,) = TT{[H,iU)~H[{t,),p{U)]ji^)} 

= J dv' (t;,oi(r',t.) -<„,(r',t,)) Tr {p(t,) [j(r),n(r')] } 

= ^Tr{p(t,)n(r)}VKoi(r,t.) -<„i(r,t,)] (A4) 
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of the system differences in time-evolution of the current expectation values. The third line 
of flA4l) results from an operator identity and by partial integration in formal manipulations 
that directly mirror those of the TD-DFT analysis. It applies because the finite (but assumed 
very large) basis for fcoi(r,t) and for ^^^[(r, t) eliminates surface contributions. 

If the potentials themselves differ at tj it follows directly that the current densities must 
differ at a time immediately thereafter. If instead the potentials only differ at some derivative 
of order A; > 1 we proceed by direct differentiation of Eq. flA4|) : 

d f .dV „ _/.5^' 



'Wtymj ^^■(r,t.) = ^Tr{p(t,)n(r)}V(^z-J [vU^.U) - v'U^.U)] ^Q. (A5) 

It follows that the current densities must differ at time f^ . 

Finally, uniqueness of the electron density n(r,t) results by direct application of the 



reductio ad absurdum argument given for infinite-system TD-DFT in Ref. 132. Using 
n(r,tj) = Tr{p(t)n(r)} and the continuity equation gives 

Qk+2 



dt^+ 



[n{v, U) - n'{v, U)] = -V • [n{r, t)Vu{r, U)] , (A6) 

gk 

u{r,ti) = —[vcoi{r, ti) -<„i(r,t,)], (A7) 



that is, in agreement with Eq. (6) of Ref. ISjJ- Use of Green's first identity shows that 

- dr u{r, U) V ■ [n(r, t)Vu] = dr n(r, t^) [Vu{r, ti)]'^ (A8) 

because our collision problem permits us to make an implicit assumption of a finite basis for 
fcoi(r, t) so that u{\r\ -^ oo) = 0. As in Refs. l3llJ32l it follows that the difference (1A6I) must be 



nonzero since n{r,ti)[Vu{r,ti)]'^ > and since n(r, tj)[V'u(r, tj)]^ cannot vanish identically. 
In summary, a reductio ad absurdum argument shows ( 01 ) that the time- dependent density 
variation n{r,t) is a unique functional of the collision potential fcoi(r, t) in the many-body 
LS collision problem that is used here to describe nonequilibrium interacting tunneling. 
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